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I.  OVERVIEW 


This  document  is  a  final  report  on  grant  #  N0001 4061 0038  (local  number  37982)  entitled  “Cooling  and 
1  rapping  of  Neutral  Atoms”  to  Harold  Metcalf  and  his  colleagues  that  terminated  30  September  2008. 

Since  our  work  on  optical  forces  on  atoms  from  non-monochroma.tic  light  began,  we  have  explored  the 
biehroniatic  force,  measured  its  properties,  provided  a  dressed  state  (Floquet)  theory,  and  developed  it  into 
a  practical  tool  for  atomic  nanofabrication.  In  addition,  we  have  explored  frequency  swept  light  in  adiabatic 
rapid  passage  to  produce  optical  forces  ten  times  larger  than  ordinary  radiative  forces  (see  Section  II  A  2). 
An  extension  of  this  theoretical  description  of  the  bichromatic  force  has  led  to  some  further  understanding 
of  the  roles  of  spontaneous  emission  and  entropy  exchange  in  all  forms  of  laser  cooling  (Section  II  B).  In 
Section  II  C  we  describe  our  successes  in  producing  our  first  nanoscale  structures.  This  is  of  special  interest 
to  the  Navy  because  of  its  possible  application  for  fabrication  of  multiply  redundant  atom  chips,  a  feature 
necessary  to  help  counteract  the  effects  of  decoherence.  Several  students  have  finished  Ph.D.’s  and  Masters 
degrees  and  there  have  been  numerous  papers  and  myriad  abstracts  for  meetings  (Section  II  D). 


II.  ACCOMPLISHMENTS  SUPPORTED  BY  THE  GRANT 
A.  Optical  Forces  in  Polychromatic  Light 

1  Introduction 

Early  laser  cooling  experiments  were  understood  in  terms  of  two-level  atoms  moving  in  a  monochromatic 
light  field  1  This  simple  view  could  describe  beam  deceleration,  dipole  force  traps,  optical  molasses,  and 
lattice's  and  band  structure  effects.  Later  the  extension  from  two-level  to  multi-level  atoms  opened  up  entirely 
new  areas  of  study  and  provided  a  whole  new  set  of  tools.  One  might  anticipate  a  similar  plethora  of  new 
phenomena  to  arise  from  the  application  of  non-monochromatic  light,  such  as  different  frequencies  or  changing 
frequencies,  but  this  topic  has  has  not  received  as  much  attention  The  first  such  experiments  were  reported 
in  Ref.  [2],  and  since  then  we  have  been  studying  optical  forces  in  such  non-monochromatic  light  [3-14). 

With  the  usual  radiative  force  from  monochromatic  light  that  has  been  exploited  for  laser  cooling  for  the 
past  25  years,  the  moment  urn  exchange  proceeds  by  absorption  followed  by  spontaneous  emission,  and  it 
saturates  at  a  maximum  value  of  Frad  =  ft/cy/2.  Moreover,  the  natural  width  imposes  a  velocity  range  for 
the  force,  limited  to  ±y//c.  In  sub-Doppler  cooling  of  multilevel  atoms,  it  is  the  dipole  force  that  works  on 
the  atoms  Even  though  the  dipole  force  does  not  saturate,  its  velocity-dependent  damping  part  is  limited 
by  7.  Also,  the  optical  {jumping  rate  7p  <  7  imposes  a  velocitv  range  ±7 p/k.  Our  measurements  have 
shown  that  polychromatic  light  doesn’t  simply  overcome  these  limits  of  force  magnitude  and  velocity  range, 
but  exceeds  them  by  an  order  of  magnitude  or  more.  In  fact,  the  limit  may  be  imposed  only  by  laser  and 
modulator  technology. 


2.  Studies  of  Adiabatic  Rapid  Passage 

Tht1  original  proposal  for  the  expiring  grant,  described  our  plans  for  exploring  and  demonstrating  extremely 
large  optical  forces  using  adiabatic  rapid  passage  (ARP),  a  long-studied  method  of  inverting  the  population 
of  a  two-level  system  that  was  well  known  since  the  early  days  of  magnetic  resonance.  Now  this  has  been 
accomplished  and  the  results  have  been  presented  in  Ref’s.  [5,  15).  We  reported  measured  forces  much  larger 
than  the  radiative  force  that  corresponded  wrell  with  these  calculations,  but  still  not  as  large  as  predicted  by 
our  model. 

Exploitation  of  ARP  for  producing  large  optical  forces  requires  repeated  frequency  sweeps  with  differently 
directed  light  beams  that  coherently  exchange  momentum  between  them,  imparting  the  difference  to  the 
atoms.  The  most  common  case  is  counterpropagating  beams  where  2 hk  is  exchanged  in  each  cycle.  The 
idea  to  investigate  ARP  grew  out  of  the  notions  of  coherent  control  of  momentum  t  ransfer  between  light  and 
atoms.  Since  the  momentum  exchange  per  cycle  2 hk  typical  atomic  momentum,  this  has  to  be  repeated 
~  1(P  times  in  order  to  have  a  significant  effect.  It  is  clear  that  ARP  can  be  effective  when  atoms  redirect 
light  from  one  beam  into  another,  and  so  the  idea  of  using  ARP  to  cause  absorption  of  light  from  one  beam 
and  stimulated  emission  into  another  came  to  mind. 
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Perhaps  the  easiest  way  to  envision  the  ARP  process 
is  in  a  dressed  atom  view  of  the  energies  of  a  two-level 
system  The  energies  of  a  pair  of  coupled  levels  in  this 
picture  are  E±  =  ±(h/2)y/S2  4-  ft2  where  6  =  u>g  -  uja  is 
the  detuning  of  the  light  at  frequency  uy  from  atomic 
resonance  at  u>a,  and  ft  ~  c(g \E  •  r\c)/h  is  the  Rabi 
frequency  that  characterizes  the  on-resonance,  electric 
dipole  interaction  between  the  light  and  atoms. 

An  important  aspect  of  the  dressed  atom  picture  for 
the  present  concern  is  the  energy  ordering  of  the  eigen¬ 
states.  In  the  low- intensity  domain  (characterized  by 
VI  <  6)  the  upper  eigenstate  approaches  the  ground  state 
| g)  and  the  lower  one  approaches  the  excited  state  |e)  for 
the  case  of  S  >  0  but  the  reverse  is  true  for  S  <  0.  A  plot 
of  E±  is  shown  in  Fig.  1  showing  these  limits  near  the 
ft  —  0  plane.  These  bare  ground  and  excited  states  |<j) 
and  |e)  are  otherwise  mixed  on  two  eigenenergy  sheets 
away  from  the  low  intensity  limit  (see  Ref.  16]). 

The*  process  of  ARP  in  this  view  involves  a  synchronized  sweep  (fast  enough  to  avoid  deleterious  effects 
from  spontaneous  emission)  of  both  the  amplitude  and  frequency  of  the  light  so  that  the  state  of  the  system 
follows  a  trajectory  similar  to  that  of  the  heavy  line  in  Fig.  1  (<$o  is  the  amplitude  of  the  frequency  sweep)  As 
long  as  travel  along  this  trajectory  is  slow  enough  to  avoid  a  non-adiabatic  transition  to  the  complementary 
eigenenergy  sheet,  the  population  will  be  completely  inverted.  Travel  on  the  lower  energy  sheet  can  be 
similarly  adiabatic  The  special  case  of  ftmax  —  <$o  was  solved  analytically  in  Ref.  [17]. 

Production  of  a  strong  optical  force  on  atoms  would  require  very  many  such  repetitive  sweeps  that  each 
exchange  momentum  hk,  and  so  the  probability  for  non-adiabatic  transitions  must  be  kept  very  small  for 
this  to  be  successful.  The  probability  of  such  unwanted  transitions  can  be  found  from  the  small  fraction  of 
population  on  the  “wrong”  energy  sheet  at  the  end  of  each  sweep  where  Vl  —  0  so  the  eigenstates  are  exactly 
the  bare  states. 


FIG.  1  A  plot  of  E±.  The  dressed  states  com¬ 
prise  two  separated  sheets  except  at  the  coni¬ 
cal  intersection  at  the  origin  The  upper  (lower) 
state  is  ground  at  ft  —  0  for  S  >  0  ((5  <  0).  The 
indicated  path  is  a  possible  trajectory  for  ARP 


FIG.  2:  A  contour  plot  of  the  numerical 
values  of  the  residual  ground  state  popula¬ 
tion  pgg  vs.  So  and  fto  after  one  sweep  (this 
is  best  viewed  in  color).  The  large  open  ar¬ 
eas  toward  the  upper  right  represent  regions 
where  pgg  <  .01  so  the  probability  of  non- 
adiabatic  transitions  is  negligibly  small.  Inter¬ 
est  lies  in  those  special  areas  that  are  close  to 
the  origin  where  pqq  is  still  tiny.  The  pathway 
along  the  vertical  axis,  =  0,  that  is  punctu¬ 
ated  by  narrow  white  regions  represents  pulses 
of  area  nn  that  also  produce  inversion  (n  — 
odd  integer).  The  indicated  points  represent 
(So/UrnMo/oJrn)  =  (a)  (2.4,  1.8),  (b)  -  (3, 

1.4),  (e)  -  (7,  7),  and  (d)  -  (14,  18).  (Figure 
from  Ref.  [3].) 


In  Ref  [3j  we  calculated  the  conditions  that  optimized  this  momentum  transfer  in  the  parameter  space  of 
So  and  ft,  by  integrating  the  equation  of  motion  of  the  Bloch  vector  on  the  Bloch  sphere.  We  very  quickly 
discovered  unexpected  behavior  of  nearly  closed  orbits  and  in  Ref.  [4]  we  provided  an  expansion  of  those 
views  to  non- ideal  sweeps.  Figure  2  shows  a  contour  plot  of  the  of  the  regions  of  parameter  space  that 
minimize  the  residual  ground  state  population  after  a  single  sweep  The  usual  viewpoint  is  that  the  ARP 
mechanism  works  in  the  large  open  areas  of  the  upper  right  quadrant,  but  there  are  clearly  other  regions 
where  it  will  also  work  well 
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We  chose*  to  do  experiments  in  the  region  near  point  (a)  of  Fig.  2  because  it  seemed  like  the  most  unusual 
place  for  ARP  to  be  effective,  and  because  it  was  the  easiest  to  reach  experimentally.  Our  successes  are 
described  in  Ref’s.  [5,  15]  and  are  summarized  in  Fig.  3.  However,  we  found  that  the  overall  magnitude  of 
the  force  was  less  than  half  of  the  value  predicted  by  our  simple  model. 


FIG.  3:  Part  (a)  shows  a  calculated  map  of  the  nonadiabatic  transition  probability  of  Ref.  [3]  from  which 
t  he  force  can  be  calculated  directly,  and  part  (b)  show's  the  force  we  extracted  from  our  measurements. 
Lighter  areas  correspond  to  larger  forces.  The  qualitative  agreement  is  quite  excellent.  (Fig.  from  Ref.  [5].) 


In  Ref.  [5]  we  suggested  that  the  problem  may  be  incomplete  optical  pumping  to  the  J  ==  1,  Mj  =  +1 
sublevel  required  to  make  this  a  true  two-level  atom  in  a  cycling  transition.  Although  the  maximum  possible 
pumping  rate  is  7/2,  our  detuning  is  sufficiently  large  that  the  excitation  rate  is  considerably  smaller.  Thus 
t  he  characteristic  time  for  even  the  first  step  in  the  optical  pumping  process  is  30  r  corresponding  to  travel 
through  nearly  1/3  of  the  interaction  region.  More  careful  modelling  of  these  processes  suggests  that  half  of 
the  atoms  are  lost,  a  result  consistent  with  our  measurements. 


We  therefore  designed  and  built  magnet  coils  and  installed  an  adequate  pre-pumping  region  upstream  of 
the  ARP  region.  The  result  was  a  substantial  increase  of  the  measured  force  and  a  paper  is  in  preparation 
[6],  but  the  magnitude  of  the  force  still  remains  well  below  that  of  our  model  Not  only  that,  but  we  found 
that  the  optical  pumping  light  beam  made  little  improvement  once  the  magnetic  environment  was  controlled. 

3.  An  Interesting  Bichromatic  Field  Result 


Many  of  t  he  features  of  the  bichromatic  force  had  already  been  explored  when  the  original  proposal  was 
written,  and  w’ere  presented  there  as  recent  accomplishments.  Almost  all  of  our  experimental  work  in  this 
area  has  previously  been  published  in  various  articles  [7  14]  so  the  topic  won’t  be  discussed  here,  However, 
we  have  found  a  number  of  curious  results  from  our  new  calculations  that  began  in  the  summer  of  2007.  One 
of  these  is  the  behavior  of  the  eigenfunctions  of  two-level  atoms  in  a  two- frequency  field.  Our  calculations  first 
reproduced  the  eigenvalues  of  the  Floquet  Hamiltonian  given  in  Ref.  Ill]  and  then  examined  the  eigenvectors 
belonging  to  them. 

Over  a  very  wide  range  of  parameters,  we  have  found  that  these  eigenvectors  conform  with  our  intuition  in 
some  cases,  but  to  give  eye-opening  results  in  others.  For  the  former,  we  note  that  the  expansion  coefficients 
have  the  expected  spatial  symmetry  and  appropriately  vanish  at.  the  nodes  of  the  fields  that  connect  them 
to  other  states.  For  the  latter,  we  were  surprised  to  see  that  the  eigenstates  of  bichromatic  light  allow'  an 
average  inverted  population  of  a  two-level  atom  in  steady  state,  a  situation  forbidden  in  monochromatic 
light .  This  surprising  result  has  now  been  described  with  a  simple  intuitive  model  on  the  Bloch  sphere  as 
well  as  formal  solution  of  the  Schrodinger  equation  in  which  the  absence  of  the  rotating  wave  approximation 
accounts  for  the  two  frequencies  [18]. 

This  result  can  be  described  in  terms  of  the  motion  of  the  Bloch  vector  R  on  the  Bloch  sphere  given  by 
dR/dt.  —  0  x  /?,  where  0  is  the  usual  artificial  3-vector  (0,0,  <S),  0  is  the  Rabi  frequency  of  the  atom-light 
transition  for  light  of  frequency  u,y  detuned  from  the  atomic  frequency  u)a  bv  6  =  -  coa.  In  such  single- 

frequency  light  the  eigenstates  (stationary  states)  are  necessarily  those  having  R  ||  ±0  so  that  dR/dt  —  0. 


These  two  dressed  eigenstates  of  the  atom-field  system  reside  on  the  upper  and  lower  hemispheres  of  the 
Bloch  sphere,  and  can  only  approach  the  equator  in  the  limit  6 /Cl  —>  0,  but  cannot  cross  it. 

Quite  the  reverse  is  true  in  a  bichromatic.  field.  The  infinite  Floquet  Hamiltonian  is  tri-diagonal  as  given 
in  Ref.  11]  and  the  spatial  dependence  of  its  eigenvalues  is  plotted  there.  However,  the  eigenvectors  were 
not  studied  until  this  past  summer  (2007)  when  a  student  plotted  their  dependence  on  Cl. 

She  found  that  in  the  limit  of  Cl/ 6  — *  0  the  eigenstates  are  the  ground  and  excited  states  as  expected, 
but.  they  become  strongly  mixed  near  Cl/ 6  ~  1,  so  we  call  them  |1)  and  |2)  in  analogy  to  the  usual  dressed 
state  labelling.  In  this  region,  the  alternate  Floquet  states  have  the  same  mixtures  of  ground  and  excited 
state  components,  again  as  expected.  However,  unlike  the  case  for  monochromatic  light  discussed  above,  her 
surprising  results  in  the  bichromatic  field  show  that  these  states  can  switch  identities  so  that  a  primarily 
ground  state  at  low  values  of  Cl/ 6  can  be  primarily  excited  near  Cl/ 6  %  2  (see  Fig.  4).  These  states  have 
stable  maxima  in  this  opposite  sense  near  Cl/6  ~  2.8. 


FIG.  4:  The  spatially  averaged  excited  state- 
component  of  type  |1)  eigenstates  vs.  Cl/ 6.  Near 
zero  it  is  mostly  ground  state,  but  becomes 
strongly  mixed  near  Cl/6  ~  1.  What  is  surpris¬ 
ing  is  the  region  above  2  where  more  than  50% 
of  its  composit  ion  is  excited  state,  which  means  it 
spends  most  of  its  time  on  the  upper  hemisphere 
of  the  Bloch  sphere. 


The  model  that  helps  us  to  describe  this  is  based  on  the  beat  frequency  pattern  of  the  bichromatic  field  as 
shown  in  Fig.  5.  We  imagine  it  to  be.  a  single  frequency  so  that  the  Bloch  sphere  picture  is  suitable,  but  its 
amplitude  is  modulated.  When  the  two  frequencies  are  equally  spaced  above  and  below  atomic  resonance, 
the  carrier  frequency  is  resonant  so  17  lies  on  the  equator,  and  its  large  amplitude  at  short  times  drives  R 
away  from  the  south  pole  very  rapidly  R  reaches  the  equator  in  a  short  time,  but  by  then  the  amplitude  of 
the  driving  field  (still  precisely  resonant)  has  begun  to  decrease.  Thus  the  traversal  of  R  across  the  upper 
hemisphere  is  much  slower,  and  in  fact,  is  a  minimum  if  it  passes  the  north  pole  just  where  the  amplitude 
vanishes  in  Fig.  5  (i.e.,  one-half  the  pattern  shown  is  a  7r-pulse.)  Then  R  spends  much  more  of  its  total  time 
in  the  upper  hemisphere  even  though  it  is  primarily  a  “ground”  state. 

FIG.  5:  The  beat  pattern  of  two  closely- 
spaced  frequencies  as  appropriate  for  the  bichro¬ 
matic  field.  It  the  early  part  of  the  beat  cycle,  the 
amplitude  is  quite  large  and  the  Bloch  vector  R 
processes  rapidly  about  17  reaching  the  equator  in 
a  short  time.  However,  its  traversal  of  the  upper 
hemisphere  is  much  slower  because  of  the  smaller 
amplitude,  so  its  dwell  time  there  is  much  longer. 


B.  Studies  of  Entropy  Exchange  in  Laser  Cooling 

1  Introduction 

Laser  cooling  is  typically  described  in  terms  of  optical  forces  from  monochromatic  light  whose  velocity 
dependence  reduces  the  width  of  the  velocity  distribution  of  an  ensemble  of  atoms.  Because  there  are  very 
many  different  phenomena  involved,  considerable  care  is  necessary  to  keep  their  descriptions  separate  from 
one  another  and  not  mix  up  the  intuitive  notions  from  each  of  these  various  effects.  For  example,  it  is  a 
long-standing  and  widely  held  tenet  in  the  laser  cooling  community  that  spontaneous  emission  is  required 
to  carry  away  the  entropy  lost  by  a  \apor  of  atoms  being  cooled.  But  this  is  often  misconstrued  to  mean 
that  spontaneous  emission  is  always  required  for  laser  cooling.  We  showed  how  cooling  could  occur  without, 
it,  and  proposed  an  experimental  test  of  this  completely  unexpected  possibility 
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We  begin  by  showing  that  spontaneous  emission  is  not  the  only  way  of  removing  the  entropy,  and  that  the 
laser  fields  themselves  are  capable  of  absorbing  it  [26].  This  is  because  spontaneous  emission  does  this  by 
redistributing  the  light  among  a  multitude  of  aecessible  states  in  the  Hilbert  space,  and  stimulated  emission 
can  do  precisely  the  same  thing.  Our  calculation  compares  the  entropy  lost  by  the  cooled  atoms  with  the 
entropy  capacity  of  the  laser  fields.  The.  description  requires  that  the  light  field  be  included  as  part  of  the 
system,  and  not  just  as  an  externally  applied  potential. 

In  the  usual  (AE,  Ap)  description  of  laser  cooling,  the  force  is  calculated  from  the  momentum  of  absorbed 
light  hvjc  =  hk  or  tr(pVEf).  Then  the  kinetic  energy  exchange  between  atoms  and  light  involves  Doppler 
shifts  or  spatially-dependent-  light  shifts  in  inhomogeneous  optical  fields.  The  velocity-dependent  force  moves 
all  the  atoms  in  an  extended  region  of  velocity  space  to  a  narrower  region  so  the  velocity  distribution  is 
compressed.  In  these  usual  (A E,  Ap)  views,  the  laser  light  is  treated  as  a  classical  field  with  a  fixed  potential, 
and  the  entropy  loss  is  usually  dismissed  with  vague  references  to  spontaneous  emission,  without  any  proof. 

When  light  is  absorbed  by  an  atom,  its  internal  energy  increase  is  compensated  by  an  energy  decrease 
of  the  light  field.  A  popular  way  to  depict  this  energy  conservation  process  is  the  Jaynes-Cummings  view, 
where  the  light  field  can  be  described  as  a  number  state.  Then  it  becomes  clear  that  some  energy  from  the 
light  field  is  transferred  to  the  atom,  but  only  if  the  light  field  is  part  of  the  system. 

2.  A  Closer  Look  at  Doppler  Cooling 

As  an  example,  consider  the  case  of  Doppler  cooling  where  the  (A  E,  Ap)  description  is  incomplete  because 
its  velocity-dependent  damping  force  does  not  conserve  energy.  The  velocity  dependence  arises  because  the 
resonance  condition,  and  hence  the  optical  force,  depends  on  the  Doppler  shift  seen  by  moving  atoms, 
u,p  =  k  ■  v,.  That  is,  the  absorbed  light  frequency  is  u'£  =  u:a  -  u:p,  where  u:a  is  the  atomic  frequency 
whereas  the  average  emitted  light  frequency  is  >  ay.  Ihe  kinetic  energy  change  on  absorption  is  AKE 
—  rnvAv  =  hkv  =  -hup  (for  v  »  Az;  =  hk/M).  Then  energy  conservation  is  satisfied  with  a  velocity- 
dependent  force  when  the  light  field  is  included  in  the  system  because  there  is  simply  an  energy  exchange 
between  the  atoms  and  the  light  through  the  Doppler  shift  ap.  Analogous  arguments  apply  to  the  velocity 
dependence  of  other  cases  of  (AE,  Ap)  exchanges. 

The  most  naive  (AE,  Ap)  view  of  laser  cooling  violates  the  unitarity  theorem.  That  is,  atoms  with  initially 
different  velocities  have  initially  orthogonal  wave  functions  because  of  their  different  deBroglie  wavelengths. 
After  cooling,  their  deBroglie  wavelengths  may  be  sufficiently  similar  that  their  wave  packets  are  no  longer 
orthogonal,  thereby  violating  the  theorem.  Unitarity  in  Doppler  cooling  may  be  rescued  by  including  the 
spontaneously  emitted  light  from  the  excited  atoms  into  the  system  because  the  fluorescence  generally 
occupies  orthogonal  states  of  the  radiation  field  thereby  preserving  the  theorem.  In  the  discussion  below,  we 
show  that  changes  to  the  laser  field  itself  are  sufficient  to  preserve  unitarity. 

In  this  review  of  Doppler  cooling,  we  have  seen  that  including  the  optical  field  provides  for  satisfying  three 
conditions:  1)  energy  conservation  between  the  atomic  internal  energy  and  the  field,  2)  energy  conservation 
between  the  atomic  motion  and  the  field,  and  3)  preservation  of  unitarity.  Similar  arguments  hold  for  other 
laser  cooling  schemes. 

3.  Phase  Space  Considerations 

We  now  extend  the  discussion  beyond  the  limits  of  the  usual  (AE,  Ap)  view  and  treat  the  laser  field  as  a 
dynamical  variable.  This  notion  is  substantiated  in  the  very  eloquent  statement  from  Ref  .  [27]: 

What  do  we  do  next?  ...  Begin  by  deciding  how  much  of  the  universe  needs  to  be  brought 
into  the  discussion.  Decide  what  normal  modes  are  needed  for  an  adequate  treatment  of  the 
problem  under  consideration.  Find  a  suitable  approximation  for  the  normal  modes;  the.  simpler, 
the  better.  Decide  how  to  model  the  light  sources  and  work  out  how  they  drive  the  wave  function 
for  the  system.’’ 

Although  this  seems  a  bit  superficial  at  first,  in  fact  it’s  very  profound.  A  complete  description  of  the 
dynamics  of  any  of  t  hese  laser  cooling  processes  requires  that  the  entire  light  field  be  considered  as  part  of 
the  system.  Only  then  can  it  absorb  entropy  and  transport  it  out  of  the  system.  Spontaneous  emission  simply 
redistributes  the  light  into  some  subset  of  a  much  larger  set  of  accessible  states,  and  stimulated  emission 
can  do  likewise.  Thus  the  exchange  of  entropy  between  the  atoms  and  the  light  field  does  not  violate  the 
Liouville  theorem,  unitarity,  or  Ref.  [28]  because  neither  the  total  entropy  of  the.  system  nor  its  phase  space 
volume  is  reduced,  but  merely  exchanged  between  its  different  parts. 


Our  1  D  comparison  with  the  entropy  capacity  of  the  laser  fields  begins  bv  first  finding  the  entropy  lost 
by  atoms,  ASa.  If  the  number  of  atoms  is  unchanged  by  the  cooling  process,  the  Sackur-Tetrode  equation 
can  be  used  to  find  A Sa  =  kg  In  where  V7^'  is  a  phase-space  volume.  (For  laser  cooling, 

AStt  <  0.)  Moreover,  the  changes  of  the  V^’s  are  expressed  by  the  product  of  the  compression  in  velocity 
space  AiJinit/Aufin*)  with  the  expansion  in  configuration  space  Axinit/Axfinai. 

For  the  latter,  we  note  that  all  laser  cooling  schemes  have  a  characteristic  cooling  length  and  time  found 
from  the  force  and  the  velocity  capture  range  vc  of  the  force.  The  natural  choice  for  the  initial  spatial  extent 
of  the  atoms  is  this  cooling  length  Ax.  The  largest  distance  atoms  can  travel  during  slowing  is  another 
Ax,  and  only  a.  few  atoms  will  reach  2 Ax.  For  calculational  convenience,  we  choose  the  final  spread  to  be 
Axfinai  %  \/2  Ax  for  each  direction,  a  total  of  2\  2  Ax  in  one  dimension.  For  the  width  of  the  initial  velocity 
distribution  Auinjt  we  take  the  velocity  capture  range  vc  and  for  Ai?finai  we  take  some  measure  of  the  cooling 
limit.  Thus  A Sa  =  kg  ln(2v'2  Aufinai/uc)  <  0,  and  is  typically  a  few  x  (-kg)  per  atom  (it’s  surprising  that 
AS a  is  so  small). 

4 .  Description  of  the  Light  Field 

It  is  straightforward  to  see  how  the  light  field  absorbs  entropy  by  calculating  its  density  matrix  pe  before 
and  after  its  interaction  with  atoms  [26,  29].  This  can  be  done  in  any  basis,  and  we  start  with  the  light 
field  in  a  pure  state  where  pe  —  {pe)2  so  the  entropy  is  zero.  For  the  interaction  with  two-level  atoms,  we 
choose  a  basis  of  number  states  for  the  light,  beams  as  in  the  Jaynes- Cummings  picture,  and  then  properly 
entangle  it  with  the  ground  and  excited  states  of  an  atom,  ]<?)  and  |e)  respectively,  to  make  a  total  wave 
function  T.  The  interaction  with  the  atoms  will  change  this  to  After  evaluating  the  matrix  p'  using 
the  operator  |'I/')('F/|  and  then  tracing  over  the  two  atomic  states  | g)  and  |e),  it  is  clear  that  p /  ^  {pe')2 
[29,  30].  Thus  the  light  field  is  no  longer  in  a  pure  state  so  its  entropy  has  increased  to  some  positive  value 
as  a  result  of  the  interaction. 

A  sample  of  atoms  immersed  in  a  light,  field  is  neither  a  closed  system  nor  is  it  in  thermal  contact  with  a 
reservoir,  so  the  ordinary  thermodynamic  entropy  cannot  be  defined.  Instead  we  use  the  information  defini¬ 
tion  S  =  kg  In  (N)  tvhere  N  is  the  number  of  states  accessible  to  the  system.  We  find  the  entropy  capacity  of 
the  light  beams  A Se  from  Ar-values  that  are  sufficiently  distinct  (small  overlap)  after  the  stimulated  emission 
processes  have  redistributed  the  light  energy  among  them. 

Although  the  natural  choice  for  a  description  of  the  laser  beams  might  seem  to  be  the  familiar  coherent 
states  |q),  the  strict  definition  of  the  |q)’s  is  not  well-suited  to  the  exchange  of  light  between  beams  caused  by 
absorption-stimulated  emission  cycles  of  atoms.  In  particular,  the  transition  term  of  the  Jaynes- Cummings 
Hamiltonian  is  oc  (a6*  +  a^6),  and  although  |a)  is  an  eigenstate  of  a,  a*|a)  is  a  complicated  object  [31-33]. 
Moreover,  the  |a)’s  are  not  eigenstates  of  the  Hamiltonian  a.^a  nor  are  they  orthogonal.  Still,  they  represent 
a  suitable  approximation  as  long  as  it’s  recognized  that  absorption  of  light  indeed  does  change  the  actual 
state  of  the  field,  even  though  in  the  exact  (ideal)  case,  |q)  is  an  eigenstate  of  a  [3 1-3-1] . 

Even  though  the  coherent  states  |q)  and  |a')  are  not  very  good  approximations,  we  may  choose  the  “distinct 
state"  criterion  from  the  overlap  formula  e~\a~oc  I  to  be  1/e.  Then  this  overlap  condition  requires  distinct 
states  to  have  n- values  that  differ  by  ±2 >Jn,  where  n  =  |q|2.  Since  Ref.  [31]  shows  that  the  approximately 
coherent  states  are  sub-Poissonnian,  there  are  actually  more  accessible  states  than  this  lower  limit  estimate. 

We  need  an  estimate  the  field  quantum  number  n,  found  from  the  amount  of  light  that  can  interact  with 
an  atom  in  a  large,  cw  laser  beam.  We  choose  a  cylinder  of  base  area  equal  to  the  on-resonance  atomic 
absorption  cross  section  a  —  3A2/2tt  and  of  length  c£cooi,  where  tcooi  is  the  the  cooling  time  mentioned 
above.  For  a  beam  of  intensity  /  =  sISAt  we  find  n  =  lotcoo\/hi'  =  s^tcoo\/2  where  /sat  =  tt/ic7/3A3.  (7 
appears  as  an  artifact  of  the  sat  uration  intensity  Js  in  /).  For  Doppler  cooling,  n  ~  103. 


5.  Entropy  Exchange 

The  maximum  change  of  n  required  to  stop  an  atom  from  v  —  vc  is  An  —  Mvc/hk  [35]  and  there  are 
many  different  values  of  An  for  the  different  atomic  velocities  in  the  sample,  so  many  states  of  the  light  field 
can  be  populated.  The  number  of  independent  states  accessible  to  the  laser  field  is  N  =  An/2\fn  so  that 
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where  vr  =  hk/M  is  the  recoil  velocity.  (Curiously,  there  is  neither  any  dependence  on  vc  in  Eq.  1  nor  buried 
in  Affinal  or  Cool-)  Moreover,  there  is  a  one-to-one  correspondence  between  the  states  of  the  light  field  and 
the  atomic  motional  states  f3G]  so  their  overlap  is  also  small.  As  long  as  77  ~  1,  meaning  A Se  ~  |ASa|,  the 
light  field  has  a  large  enough  capacity  to  absorb  the  entropy  lost  by  the  atoms. 

Further  evaluation  and  discussion  of  Eq.  1  must  be  done  case  by  case  and  the  results  are  summarized  in 
Table  1  [26].  The  first  two  entries  for  Doppler  molasses  are  readily  calculated  from  the  well-known  Doppler 
limit  and  velocity  damping  constant  [1,  371.  The  cooling  limit  for  Sisyphus  cooling  by  polarization  gradients, 
typically  a  few  x  vr  1  38],  is  taken  from  Eq.  4.37  of  Ref.  [39]  ((5  =  -  u;a).  The  second  entry  is  readily 

calculated  from  the  increased  velocity  damping  coefficient  of  Sisyphus  cooling  [39,  40]. 

The  first,  entry  for  the  bichroma  tic  force,  whose  two  frequencies  are  detuned  by  ±  <5,  comes  from  as¬ 
suming  that  the  atoms  are  finally  distributed  between  two  ground  states  (it  could  be  more)  [11],  and 
the  second  entry  comes  from  dividing  the  velocity  cooling  range  by  the  (approximately  constant)  accel¬ 
eration  2hk6/nM  [41].  The  value  s  =  3$2/72  optimizes  the  bichromatic  force  (typically  6  ~  40">). 


TABLE  I.  A  summary  of  the  values  for 
entry  into  Eq.  1.  The  velocity  limit 
of  Doppler  cooling  is  v&  =  \fhry/(2M). 
See  text,  for  a  discussion.  For  Sisyphus 
cooling  S  =  lu'(  —  u>q  is  the  laser  detuning. 
For  the  bichromatic  case,  2 S  is  the  fre 
quency  difference  between  the  two  fields, 
and  typically  6  ~  40y. 

Clearly  77  ~  1  for  all  three  cases  so  we  can  conclude  that  the  laser  beams  themselves  have  sufficient  capacity 
to  absorb  the  entropy  lost  bv  the  cooled  atoms  by  redistribution  of  the  laser  light  into  a  larger  number  of 
accessible  states.  Since  the  system  is  uopen’\  the  outgoing  light  beams  carry  away  the  entropy.  Spontaneous 
emission  is  not  required  for  this  aspect  of  laser  cooling.  In  some  sense,  this  may  be  related  to  cavity  cooling 
(see  Ref.  [42]  and  references  therein). 

We  see  that  the  entropy  lost  by  the  atoms  can  be  transferred  to  the  light  field,  and  it’s  not  dissipated  until 
the  outgoing  light  hits  the  walls.  Since  the  walls  are  NOT  part  of  the  system,  this  final  destruction  of  the 
light  field  is  indeed  a  dissipat  ive,  non-conservative,  and  irreversible  process.  The  walls  are  not  part  of  the 
system,  just  as  the  empty  modes  into  which  spontaneous  emission  dumps  the  light  are  not  part  of  the  system 
in  the  common  ( A 7T,  Ap)  pictures  of  Doppler  cooling.  The  model  here  is  different  from  those  of  previously 
studied  cases  [28,  43.  44]  because  the  applied  light  field  is  a  dynamical  variable  of  the  Hamiltonian. 
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6.  The  Hole  of  Spontaneous  Emission  in  Laser  Cooling 

Although  it’s  now  clear  that  spontaneous  emission  is  not  required  to  carry  away  the  entropy  removed  from 
a  vapor  atoms  during  laser  cooling,  it  is  also  true  that  spontaneous  emission  is  indeed  crucial  for  energy 
exchange  in  some  forms  of  laser  cooling.  Since  excited  atoms  don’t  absorb  resonantly,  their  return  to  the 
ground  state  without  the  spontaneous  emission  needed  for  further  absorption  can  only  occur  by  stimulated 
emission.  If  it  occurs  by  the  original  exciting  beam,  there  is  no  net  energy  exchange  in  the  process.  If 
stimulated  emission  is  induced  by  another  beam,  for  example  the  counterpropagating  beam  in  Doppler 
molasses  or  Sisyphus  cooling,  then  the  absorption-stimulation  emission  sequence  in  the  opposite  order  is 
equally  likely  when  spatially  averaged  over  a  wavelength,  and  on  average  there  is  no  net  energy  transfer 

For  Doppler  molasses,  only  spontaneous  emission  between  absorptions  can  exploit  the  Doppler  shifts  to 
allow  energy  loss  because  the  average  frequency  of  the  emitted  light  is  ljq  >  aj?.  For  Sisyphus  cooling, 
spontaneous  emission  is  needed  for  the  t  ransitions  between  the  top  of  one  potential  hill  of  the  light  shift  to 
the  bottom  of  another  [39].  Without  it,  the  energy  exchange  is  limited  to  the  light  shift,  and  that  is  usually 
much  smaller  than  atomic  kinetic  energies.  Thus  the  indispensable  role  of  spontaneous  emission  in  many 
forms  of  laser  cooling  is  not  for  entropy  dissipation,  but.  for  energy  transfer. 

By  contrast.,  the  bichromatic  force  (BF)  has  no  such  limits  [2,  11].  With  appropriate  parameters  and 
spatial  offsets  of  the  standing  waves  of  the  two  different  frequencies,  there  are  positions  where  the  light  shifts 
and  detuning  just  cancel,  and  exact  crossings  of  the  eigenstates  can  occur.  Atoms  can  mediate  the  exchange 
of  red- detuned  light  for  blue,  and  vice  versa,  resulting  in  energy  exchange  between  the  atoms  and  the  light 
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beams  1 1  Thus  energy  conservation  and  entropy  dissipation  for  the  BF  are  accomplished  differently  so 
that  spontaneous  emission  is  not  necessary  for  any  of  them.  This  mechanism  lias  no  significant  limit  on  the 
energy  scale  of  interest  here,  and  more  important  has  no  counterpart,  in  monochromatic  light  26]. 

The  BF  is  indeed  velocity  dependent,  hut.  it’s  not  at  all  like  the  more  familiar  friction  forces  of  other 
laser  cooling  schemes.  Here  the  non- adiabatic  transitions  provide  the  origin  of  the  velocity  dependence  of 
the  force,  and  depend  exponentially  on  the  velocity  |11].  There  is  no  traditional  friction  or  damping  force. 
Thus  (he  BF  is  different  from  other  laser  cooling  schemes  because  spontaneous  emission  is  NOT  needed  to 
mediate  the  energy  exchange  nor  to  produce  a  velocity-dependent,  force. 

We  have  seen  how  the  light  beams  have  the  capacity  to  absorb  the  entropy  lost  by  laser  cooled  atoms, 
and  that  spontaneous  emission  is  not  required  for  this  task.  We  have,  also  seen  that  it  is  indeed  required  to 
mediate,  the  energy  loss  of  the  atoms  in  some  forms  of  laser  cooling.  The  combination  of  these  two  results 
provides  a  completely  different  view  of  the  role  of  spontaneous  emission  in  laser  cooling. 

C.  Nanofabrication  with  He*  Using  Bichromatic  Beam  Collimation 

1.  Introduction 

Atomic  nano-fabrication  (ANF)  has  many  special  features  that  distinguish  it  from  other  methods  [19] 
and  these  were  outlined  in  the  proposal  for  the.  expired  grant.  At  that  time  we  proposed  to  do  ANF  with 
metastable  23S  helium  (He*)  very  efficiently  because  we  could  use  the  bichromatic  force  to  get  a  high  He* 
intensity  so  that  the  exposure  time  could  be  very  short.  Here  we  report  success  in  doing  precisely  this. 

The  ability  to  record  and  preserve  the  positions  of  atoms  in  a  permanent  structure  on  a  sub-micron  scale 
has  myriad  possible  applications.  To  begin,  the  structure  itself  may  be  of  value.  Such  nano-scale  patterns 
can  be  made  with  precise  long-range  regularity  whose  spacing  is  known  to  spectroscopic  precision.  Examples 
include  gratings  for  short  wavelength  radiation,  photonic  crystals,  metamaterial  samples,  and  many  others. 
Second,  the  distribution  of  atoms  is  determined  by  light  fields  so  that  the  pattern  produced  can  be  used 
as  a  measure  of  the  light  intensity  distribution.  For  example,  very  small  deviations  caused  by  wave  front 
imperfections  could  be  detected.  Finally,  the  process  by  which  the  light  fields  influence  the  atomic  motion 
in  forming  such  patterns  is  an  interesting  topic  of  study  by  itself. 

Atom  lithography  with  a  positive  resist  works  by  destroying  the  chemical  bonds  in  a  layer  of  polymeric 
molecules  of  self-assembled  monolayers  (SAM’s),  then  dissolving  the  damaged  molecules,  and  finally  etching 
the  exposed  area  (see  Fig.  6).  Using  the  metastable  23S  state  of  helium  (He*)  is  ideal  for  this  because  its 
internal  energy  of  20  eV  is  higher  than  any  other  metastable  atom.  Thus  it  is  most  effective  for  exposing 
a  resist  with  minimum  dosage  and  therefore  shortest  exposure  time  and  thus  minimum  restrictions  on  the 
atomic  beams  and  laser  control  [19]. 


FIG.  6:  The  drawing  shows  an  artistic  rendition  of 
a  fabricated  SAM  on  a  substrate.  The  long  chain 
molecules  orient  themselves  with  their  hydropho¬ 
bic  heads  bound  to  the  substrate  and  their  long- 
chain  hydrophillic  tails  sticking  out  into  the  solu¬ 
tion.  The  result  resembles  an  ordinary  carpet  pile  as 
shown.  Incident  He*  atoms  release  theii  20  eV  inter¬ 
nal  energy  somewhere  along  the  chain  and  therefore 
weaken  or  actually  break  the  bond  that  holds  it  to¬ 
gether.  Then  an  acid  solution  dissolves  away  the 
alcanet  hides  and  leaving  the  substrate  exposed  for 
further  etching. 


In  our  first  experiments,  we  planned  to  pattern  a  resist  on  a  SAM  over  gold  using  a  physical  mask  (just 
a  screen)  with  our  bright  He*  beam.  We  found  excellent  results  with  a  beam  that  was  collimated  with 
the  bichromatic  force  followed  by  two  transverse  Doppler  molasses  stages  for  further  cooling.  The  metal 
screen  projected  its  image  on  a  SAM  of  nonanethiol  and  the  undisturbed  SAM  regions  then  protected  a  gold 
coated  silicon  wafer  from  a  wet  chemical  etch.  The  open  areas  of  the  screen  allowed  incident  He*  to  damage 
the  SAM  molecules  by  depositing  their  20  eV  of  internal  energy  on  the  surface.  Samples  created  with  this 
method  had  an  edge  resolution  of  about  60  inn  that  was  measured  with  an  atomic  force  microscope. 
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Then  we  succeeded  with  nanolithographv  using  a  standing  wave  to  make  the  patterns.  Moreover,  the 
parameters  of  this  optical  standing  wave  extended  into  a  previous  unexplored  range.  In  particular,  there  are 
two  domains  of  laser  intensity  used  for  steering  the  atoms  into  the  desired  patterns  Our  experiments  are  in 
the  lower  intensity  region  called  the  focusing  regime  (at  higher  intensities  it’s  called  the  channelling  regime). 
Although  our  numerical  simulations  of  atomic  focussing  by  a  laser  standing  wave  suggest  that  we  cannot 
produce  patterns  in  this  regime  when  the  atomic  beam  has  a  relatively  large  angular  divergence,  we  have 
collimated  our  beam  sufficiently  to  accomplish  such  patterning  at  such  lower  intensities.  Smaller  features 
can  be  obtained  in  the  focussing  regime  than  in  the  channelling  regime  [20]. 

2  Numerical  Calculations 

In  preparation  for  our  experiments,  we  did  numerical  calculations  to  determine  the  neighborhood  of  our 
operating  parameters  The  first  of  these  were  simple  trajectories  of  atoms  incident  perpendicular  to  a 
standing  wave  with  transverse  Gaussian  intensity  distribution.  We  found  trajectories  corresponding  to  both 
the  focusing  and  channelling  regimes,  and  especially  had  a  close  look  at  the  region  between  them. 

Since  our  choice  of  a  perfectly  collimated,  mono- velocity  beam  incident  exactly  at  7t/2  is  unrealistic,  we 
also  performed  a  Monte-Carlo  simulation  of  the  paths  and  landing  points  of  the  atoms.  The  single  longitudi¬ 
nal  and  transverse  velocities  were  replaced  by  distributions  characterized  by  widths  that  are  correlated.  The 
simulation  results  showed  the  formation  of  half- wavelength  size  structures  over  a  wide  range  of  parameters. 
I 'he  pure  focusing  regime  for  the  parameters  reflecting  our  experimental  setup  was  not  found  to  produce 
patterning;  but.  the  simulations  indicated  that  pattern  formation  was  more  sensitive  to  atomic  beam  diver¬ 
gence  as  opposed  to  atomic  beam  chromaticity  In  both  the  trajectory  and  Monte-Carlo  calculations,  our 
results  differed  very  little  from  those  previously  published  in  Ref’s  [21,  22]. 

3.  Experiment 

We  achieved  neutral  atom  lithography  using  a  bright  beam  of  He*  that  originates  from  a  reverse  flow  DC 
discharge  source  [23,  24].  It  has  a  slightly  supersonic  longitudinal  velocity  distribution  centered  near  1100 
m/s  and  is  about  ±  200  m/s  wide.  This  atomic  beam  is  collimated  with  the  bichromatic  force  [12],  followed 
by  three  optical  molasses  velocity  compression  stages  [25].  Because  bichromatic  collimation  makes  such  an 
intense  lie*  beam,  our  exposure  time  is  measured  in  minutes  instead  of  hours.  We  have  observed  the  focusing 
and  channeling  of  t  he  He*  beam  by  the  dipole  force  the  atoms  experience  while  traversing  a  standing  wave 
of  A  —  1083  nrn  light  tuned  490  MHz  above  the  23S !  — ►  23?2  transition. 

In  our  experiment  the  bichromatic  force  is  implemented  with  counter-propagating  light  beams,  each  con¬ 
taining  two  frequencies  that  are.  separated  by  26  =  2tt  x  120  MHz.  All  the  light  originates  from  diode 
lasers  stabilized  by  saturated  absorption  spectroscopy,  frequency  shifted  by  various  AO  M’s,  and  amplified 
by  Yb-doped  fiber  amplifiers  The  laser  linewidths,  both  before  and  after  the  fiber  amplifiers,  are  measured 
by  heterodyne  spectroscopy  to  be  less  than  0.25  MHz 

If  frequencies  of  these  beams  were  centered  about  the  atomic  resonance,  the  dependence  of  the  force  on 
atomic  velocity  v  would  be  symmetric  about  v  —  0  over  the  velocity  range  ±6 /2k  and  nearly  vanish  elsewhere 
[2]  Since  collimation  requires  the  force  to  be  antisymmetric  about  v  =  0,  we  shift  the  velocity  dependence 
by  6/2k  without  changing  its  shape  by  shifting  the  laser  frequencies  appropriately  [12].  The  resulting  force 
is  unidirectional  so  there  are  two  collimation  regions  for  each  of  two  dimensions,  making  four  bichromatic 
force  regions.  The  four  sequential  bichromatic  force  regions  are  placed  as  close  as  possible  to  the  cone  of 
He*  atoms  emerging  from  the  source  so  that  the  collimated  beam  has  a  minimum  diameter  [25].  Each  region 
is  11  mm  long,  and  their  apodized  Gaussian  beam  profiles  carry  an  average  intensity  of  4000  /s,  where  the 
saturation  intensity  Is  =  tt/ic/3A3t  %  0.16  mW/cm2. 

The  transverse  velocity  spread  of  the  atoms  in  the  beam  after  the  bichromatic  force  regions  is  well  above 
±10  rn/s  corresponding  to  ~  ±9  mrad,  and  is  not  suitable  for  our  purposes.  Therefore  we  have  a  ‘‘booster” 
molasses  stage  with  detuning  many  times  larger  than  the  natural  width  7  to  capture  atoms  in  this  large 
angular  spread,  followed  by  an  ordinary  Doppler  molasses  to  bring  the  atomic  beam  divergence  down  to 
nearly  1  mrad  with  very  little  loss  of  atoms  (the  Doppler  limit  is  ~  0  25  mrad).  This  collimation  occupies  8.2 
cm  of  our  beamline  and  delivers  ~  1.5  x  109  atoms/s-mm2  at  our  sample,  68  cm  from  t.he  He*  source  [25]. 

Our  beam  of  He*  atoms  is  focussed  into  lines  by  a  standing  wave  field  of  A  —  1083  nm  light  tuned  ~ 
490  MIIz  above  atomic  resonance  (atoms  attracted  toward  the  nodes).  It  has  an  elliptical  cross  section  with 
waist  of  ~  15  mm  parallel  to  the  substrate  surface  and  of  ~  330  fim  along  the  atomic  beam  path  and  a  peak 
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intensity  of  1.48  W/em2.  The  intensity  of  this  standing  wave  can  be  varied  from  the  intensity  of  the  focusing 
regime  to  over  37  W/em2,  well  into  the  channelling  regime  [22]. 

Our  samples  are  built  oil  single-crystal  silicon  wafers  that  have  a  200  A  layer  of  gold  evaporated  onto 
their  1 00]  surface  over  a  5  A  chromium  adhesion  layer.  These  wafers  are  grown  and  coated  commercially. 
The  SAM  is  assembled  onto  the  gold  surface  by  submersing  the  wafer  in  a  1  mM  solution  of  nonanethiol  in 
ethanol  for  13  to  20  hours.  The  long  chain  molecules  orient  themselves  with  their  hydrophobic  heads  bound 
to  the  gold  substrate  and  their  long-chain  hydrophillic  tails  sticking  out  into  the  solution  as  in  Fig.  6. 

The  incident  He*  attacks  the  structure  of  the  SAM  molecules  thereby  making  it  soluble  in  a  wet  chemical 
etch  where  the  undisturbed  SAM  protects  parts  of  the  gold  layer  and  the  unprotected  gold  regions  are 
dissolved.  Thus  a  10  minute  exposure  deposits  one  He*  in  each  ~  100  A2  area  without  any  focussing,  and 
the  focussing  could  easily  increase  this  dosage  to  a  few  He*  per  SAM  molecule  site  (area  ~  30  A2,  estimated). 
Samples  created  with  this  method  have  an  edge  resolution  of  ~  80  nm  that  we  measured  with  an  atomic  force 
microscope.  This  seems  to  be  limited  by  the  domain  granularity  of  the  gold  layer  and  the  etching  process. 
This  method  of  neutral  atom  lithography  is  analogous  to  ordinary  resist-based  technologies  that  are  used  in 
most  conventional  lithography  processes. 

Figure  7  shows  ~  200  lines  spanning  about  100  x  65  /xm  of  a  3  x  3  mm  region  (less  than  1/1000  of  the 
total  area)  that  was  quite  uniformly  patterned.  The  unfocused  dosage  of  He*  was  3x  1012atoms/mm2  during 
the  36  minute  exposure.  The  intensity  of  the  light  mask  beam  used  for  patterning  was  5.93  W/cm2,  about 
four  times  /„  and  corresponding  to  ~  4  times  the  1.48  W/cm2  threshold  for  the  focusing  regime.  It  is  cut 
in  the  middle  by  the  substrate,  as  on  the  left  side  of  Fig.  3  in  Ref.  [19].  Thus  its  width  is  3  mm  along  the 
surface  perpendicular  to  its  /c-vectors  (twice  the  waist)  and  the  atoms  travel  ~  330  /im  through  half  of  it 
before  hitting  the  SAM-coated  substrate. 


FIC.  7:  Scanning  electron  microscope  image 
of  one  sample.  There  are  ~  200  vertical  lines  on 
this  figure  separated  by  541.3  nm  =  A/2.  The 
unfocused  peak  dosage  was  3xl012  atoms/inm2 
during  the  31  minute  exposure.  The  inset  is 
a  small  portion  magnified  by  5  since  aliasing 
of  the  image  in  print  or  on  screen  may  cause 
deception  even  though  its  texture  has  been 
smoothed  for  reproduction. 
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FIG.  8:  Atomic  force  microscope  image  of 
the  same  wafer  as  Fig.  7.  This  image  has  a  5 
pt  smoothing  but  still  clearly  shows  both  the 
graininess  of  the  original  gold  coating  and  the 
resolution  of  the  AFM. 


Figure  <8  shows  the  AFM  image  of  the  same  sample.  The  angular  view  shows  that  the  depth  of  the  lines 
is  nearly  the  full  10  nm  thickness  of  the  gold  layer  and  that  the  spatial  resolution  is  limited  by  its  graininess 
to  about  1/5  of  the  line  separation,  100  nm. 

Our  fabrications  have  been  reliable  and  repeatable,  sample  after  sample,  over  a  period  of  months.  In  a 
standing  wave  arranged  to  focus  the  He*  into  lines  separated  by  A/2  on  the  sample,  our  lines  cover  the  entire 
exposed  length  of  the  substrate,  about  3  inm.  They  are  3  mm  long,  corresponding  to  about  twice  the  beam 
waist  of  the  laser  standing  wave.  Thus  there  are  ~6  x  103  lines  of  length  ~  2800  A  as  shown  in  Fig.  7  [25]. 
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D.  Student  Progress  and  Other  Contributions  of  Naval  Interest 

Since  our  institution  is  a  university,  our  primary  outcome  of  special  interest  to  the  Navy  is  educational. 
Our  program  cannot  possibly  compete  with  major  institutions  devoted  to  time- keeping  and  navigation,  but 
we  DO  provide  the  educated  scientists  necessary  for  their  operation.  For  example,  current  graduate  student 
Matthew  Eardley  is  pursuing  his  Ph.D.  thesis  research  in  the  Chip  Scale  Atomic  Devices  group  at  N.I.S.T 
in  Boulder,  CO.  The  group  is  developing  gyroscopes,  magnetometers,  and  clocks  that  all  take  advantage  of 
MEMS  technology  to  produce  low-cost,  low'-pow'er,  portable  microchip-scale  devices  that  derive  high  stability 
and  accuracy  from  the  energy  levels  of  laser-pumped  atoms  in  microfabricated  vapor  cells.  Although  this 
kind  of  w'ork  cannot  be  supported  in  our  laboratory,  the  education  and  training  he  received  here  is  vital  to 
his  thesis  research  there. 

As  another  example.  2002  Ph.D.  Matt  Cashen  did  a  post-doc  with  Mark  Kasevich  at  Stanford  University 
on  precision  gravimetry,  certainly  important  for  submarine  detection.  Then  he  started,  a  permanent  job  in 
the  Lasers  Division  of  Raytheon  Corp.  wrorking  on  atomic  clocks.  The  results  of  his  research  has  been  passed 
ofT  to  Syrnrnetricom  for  further  development.  Matt  has  used  his  background  to  advance  the  state  of  the  art  in 
sensors  for  precision  inertial  navigation  systems  such  as  accelerometers  and  gyroscopes.  All  such  devices  are 
based  on  ultra.-cold  atomic  techniques  and  have  the  potential  to  revolutionize  military  navigation  systems. 
Scientists  with  his  kind  of  training  in  laser  cooling  and  trapping  techniques  are  vital  to  the  development  of 
next  generation  sensors  for  military  applications.  This  shows  the  value  of  training  Ph.D.  scientists  who  will 
enable  the  transition  of  this  technology  from  the  lab  into  the  field. 

There  are  other  facets  of  our  educational  activities,  and  part  of  its  graduate  component  is  summarized 
in  the  chart  below.  In  the  past  three  years  we  will  have  finished  four  Ph.D.’s,  three  of  whom  have  already 
taken  jobs  in  physics.  The  six  visiting  M.A.  students  from  our  exchange  programs  (five  from  Wurzburg) 
have  all  returned  home  and  begun  Ph.D.  studies  at  various  universities.  Some  of  our  undergraduates  have 
won  prestigious  awards,  and  all  of  them  are  in  Ph.D  programs  at.  excellent  universities.  Our  Laser  Teaching 
Center  students  are  not  listed  in  this  chart 


Year 

Ph.D. 

M.A. 

L^ndergraduate 

2004 

Matt  Part  low 
post-doc  at  U.  Toronto 

Jorg  Bochrnann 

Univ.  of  Wurzburg,  now  at  Miinchen 

Melissa  Friedman1 
now  grad  at  Oxford 

2005 

Matthias  Riedmann 

Univ.  of  Wiirzburg,  now  at.  Hannover 

Kyung  Choi 

Stony  Brook  under  grad 

2006 

Seung-Hyun  Lee 

now  at  Korean  Military  Institute 
Xiyue  Miao 

postdoc  at  Stony  Brook 

Esther  Wertz 

returned  to  Uruv.  Parvs  VIJ 

Mike  Keller 

Univ.  of  Wiirzburg ,  now  at  Vienna 
Andy  Vernaleken 

Univ.  of  Wurzburg,  now  at  Miinchen 

Maaneli  Derekhshani 
Stony  Brook  under  grad 
Kvung  Choi 
now  grad  at  Caltech 
Tak  Chu  Li 
now  grad  at  U.  MD. 

2007 

Benedikt  Scharfenberger 

Univ.  of  Wiirzburg ,  now  at  Sitgcn 

Thien  An  Nguyen 
Stony  Brook  under  grad 

Won  prestigious  Marshall  Fellowship  and  now  in  graduate  school  at  Oxford  University. 


1 


13 


* 


[1]  H.  Metcalf  and  P.  van  der  Strat.en,  Laser  Cooling  and  Trapping  (Springer,  New  York.  1999). 

[2]  ,L  Soding,  R.  Griinin,  5  Ovchinnikov,  P.  Bouyer,  and  C.  Salomon,  Phys.  Rev.  Lett.  78,  1420  (1997). 

[3]  T.  Lu  X.  Miao,  and  H.  Metcalf,  Phys.  Rev  A  71  R061405  (2005) 

[4]  T.  Lu,  X.  Miao,  and  H.  Metcalf,  Phys.  Rev  A  75,  063422  (2007). 

[5]  X.  Miao,  E.  Wertz,  M.  G.  Cohen,  and  H  Metcalf,  Phys.  Rev.  A  75  011402  (2007). 

[6]  X  Miao,  B.  Scharfenberger,  M.  G.  Cohen,  and  H.  Metcalf  (????),  in  preparation. 

[7]  M.  Williams,  F.  Chi,  M.  Cashen,  and  H.  Metcalf,  Phys.  Rev.  A  60.  R1763  (1999). 

[8]  M.  Williams,  F.  Chi,  M.  Cashen,  and  H.  Metcalf,  Phys.  Rev.  A  61,  023408  (2000). 

[9]  M.  Cashen  and  H.  Metcalf,  Phys.  Rev,  A  63,  025406  (2001). 

[10]  M.  Cashen,  O.  Rivoire,  V.  Romanenko,  L.  Yatsenko,  and  H.  Metcalf,  Phys  Rev.  A  6  4,  063411  (2001). 

[11]  L.  Yatsenko  and  H.  Metcalf,  Phys.  Rev.  A  70,  063402  (2004). 

[12]  M.  Partlow,  X.  Miao  J.  Bochinann,  M.  Cashen,  and  H.  Metcalf,  Phys.  Rev.  Lett.  93.  213004  (2004). 

[13]  M.  Cashen,  O.  Rivoire,  L.  Yatsenko,  and  H.  Metcalf,  J.  Opt.  B:  Quant.  Semiclass.  4,  75  (2002). 

[14]  M.  Cashen  and  H.  Metcalf,  J.  Opt.  Soc  Am.  B  20,  915  (2003). 

[15]  X.  Miao,  Ph.D.  thesis,  Stony  Brook  University  (2006). 

[16]  L.  P.  Yatsenko,  S.  Guerin,  and  H.  R.  Jauslin,  Phys.  Rev.  A  65,  043407  (2002). 

[17]  D.  Sawicki  and  J.  H.  Eberly,  Optics  Express  4,  217  (1999). 

[18]  Paul  Bennan,  private  comiri  lnication. 

[19]  D.  Meschede  and  H.  Metcalf,  J.  Phys.  D:  Appl.  Phys.  36  R17  (2003). 

[20]  W.  R  Anderson,  C.  C.  Bradley,  J.  J.  McClelland,  and  R.  J.  Celotta,  Phys.  Rev.  A  59,  2476  (1999). 

[21]  J.  J  McClelland,  R.  E.  Scholten,  E.  C.  Palm,  and  R.  Celotta,  Science  262,  877  (1993). 

[22]  J.  J.  McClelland,  J.  Opt.  Soc  Ain.  B  12,  1761  (1995). 

[23]  J.  Kawanaka,  M.  Hagiuda,  K.  Shimizu,  F.  Shimizu,  and  H.  Takurna,  Appl.  Phys.  B.  56,  21  (1993). 

[24]  H  Mastwijk,  M.  van  Rijnbach  J.  Thomsen  P.  van  der  Straten,  and  A.  Niehaus,  Eur.  Phys.  J.  D  4,  131  (1998). 

[25]  C.  Allred,  Ph.D.  thesis,  Stony  Brook  University  (2009). 

[26]  H  Metcalf,  Phys.  Rev.  A  77,  R  061401  (2008)" 

[27]  W  E  Lamb,  App.  Phys.  B  60,  77  (1995). 

[28]  W.  Ketterle  and  D.  E.  Pritchard,  Phys.  Rev.  A  46,  4051  (1992). 

[29]  E.  Boukobza  and  D.  J.  Tannor,  Phys.  Rev.  A  71,  063821  (2005). 

[30]  D.  Tannor  and  A.  Bartana,  Journal  of  Physical  Chemistry  A  103,  10359  (1999). 

[31]  G.  S.  Agarwal  and  K.  Tara,  Phys.  Rev.  A  43,  492  (1991). 

[32]  A  Zavatta,  S.  Viciani,  and  M.  Bellini,  Science  306,  660  (2004). 

[33]  Y.  Parigi,  A.  Zavatta,  M.  Kiin,  and  M.  Bellini,  Science  317,  1890  (2007) 

[34]  B.  R.  Mollow,  Phys.  Rev.  A  12,  1919  (1975). 

[35]  The  momentum  change  for  an  absorption-stimulated  emission  cycle  between  counterpropagating  beams  is  2/i/c, 
but  two  beams  are  affected  so  the  factors  of  2  cancel. 

[36]  S.  J.  van  Enk  and  H.  J.  Kimble,  QIC  2.  1  (2002),  URL  http:  //vwv.  citebase  .  org/abstract7ids=oai :  arXiv  .  org: 
quant -ph/01*/, 07088. 

[37]  D.  Wineland  and  W .  Itano,  Phys.  Rev.  A  20,  1521  (1979). 

[38]  C.  Salomon,  J.  Dalibard  W.  Phillips,  A.  Clairon,  and  S.  Guellati,  Europhvs.  Lett.  12  683  (1990). 

[39]  J.  Dalibard  and  C.  Cohen-Tannoudji.  J.  Opt.  Soc.  Am.  B  6  2023  (1989). 

[40]  Caution  is  needed  since  Eq.  4.3  of  Ref.  [39]  defines  the  Rabi  frequency  to  be  twice  as  large  as  that  used  here. 
Our  definition  is  Mi  =  d  •  £. 

[41]  In  Ref’s.  [9,  12]  the  velocity  range  was  measured  to  be  S/k  (±6/2k)  but  in  Ref.  [11]  it  was  shown  that  the  part, 
of  the  bichromatic  force  that  arises  from  the  mechanism  described  therein  has  a  range  S/2k  (±6/4k). 

[42]  K.  Mfirr,  Phys.  Rev  Lett.  96  253001  (2006). 

43]  M  E.  Carrera- Patio  and  R.  S.  Berry,  Phys.  Rev.  A  34,  4728  (1986). 

44]  S.  van  Enk  and  CL  Nienhuis,  Phys.  Rev.  A  46,  1438  (1992). 


